Electrostatic pair creation and recombination in quantum plasmas 
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The collective production of electron-positron pairs by electrostatic waves in quantum plasmas is 
investigated. In particular, a semi-classical governing set of equation for a self-consistent treatment 
of pair creation by the Schwinger mechanism in a quantum plasma is derived. 



I. INTRODUCTION 
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In quantum field theory, and quantum electrodynam- 
ics (QED) in particular, the vacuum is no longer trivial, 
but may under certain circumstances act as a medium. 
Effects such as elastic light-by-light scattering as pho- 
ton splitting are well known examples of the influence 
of the nonlinear quantum vacuum The concept of 
pair creation is of great interest both from a principal 
point of view, as well as for near-future applications. 
In particular, the production of pairs in an intrinsically 
nonperturbative process, and as such poses new require- 
ments for its theoretical treatment. Lately, there has 
been much interest in pair creation due to nontrivial gen- 
eralizations of the constant electric field case (see e.g. 
HHHIlllllllli and references therein), as con- 
sidered by Schwinger for quantum electrodynamics |lOj |. 
The next generation laser systems, such as the X-ray free 
electron laser, will be able to produce pairs, and this ne- 
cessitates the study of collective and self-consistent pair 
production effects. Thus, in this spirit, the effects of pair 
production on circularly polarized light was considered in 
Refs. ^3 and 0, and it was found that electromagnetic 
degrees of freedom were dissipated into pairs in a self- 
consistent manner. However, it is not unreasonable to 
assume that within laser-plasma and astrophysical sys- 
tems, strong electrostatic fields may be built up, thus 
posing a new problem concerning pair creation. 

In this paper we will study the effects of electrostatic 
pair creation in a three component plasma, consisting 
of ions, electrons, and positrons. The quantum proper- 
ties of electrons and positrons are partially included by 
using a semi-classical approximation of the many-body 
Schrodinger system. The effect of the latter is to intro- 
duce higher order dispersion in the electron/positron mo- 
mentum equations. The combined effect of pair creation 
and quantum plasmas can thus be studied using and ef- 
fective semi-classical theory in terms of macroscopic fluid 
variables. 
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Here p is the electron/positron microscopic momentum, 
E is the electric field strength, e is the magnitude of 



Ic^/eU - IQi^V cm- 



is 



the electron charge, i^crit = 
the Schwinger critical field, A = h/mgC is the Compton 
wavelength, rUe is the electron rest mass, c is the speed 
of light in vacuum, U in Planck's constant divided by 27r, 
and the momentum distribution function is normalized 
according to / Fd^p = 1. produced at rest, thus the 
delta function in momentum space, and neglected the 
small correction due to the electron/positron momentum 
orthogonal to the electric field. We note that the expres- 
sion 1^ should be used with caution, as pointed out in 
Ref. j^], as the electrons/positrons does not experience 
the local field. This can be partially remedied by averag- 
ing the field over a suitable spatial volume. Here we will 
at present consider the local field approximation, keep- 
ing in mind the above. Apart from the pair creation rate, 
there will also be a pair recombination, or annihilation, 
rate given by = aVcUp [vp = crvcne) for the electrons 
(positrons), where a is the cross section, Vc is a charac- 
teristic velocity of the recombination, and (rip) is the 
electron (positron) number density. 

For electrostatic oscillations, the electric field is deter- 
mined by Poisson's equation 



V ■ E = {-eue + erip + ZiCUi + Ppoi)/eo- 



(2a) 



where (up) is the electron (positron) number density 
in the laboratory frame, Ui is the ion number density, ZiC 
is the ion charge density and ppoi is a vacuum polarization 
charge density. Ampere's law gives 
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II. THEORY 



For fields varying slowly compared to the Compton 



frequency We 



' /fi, where me is the electron rest 



mass, c is the velocity of light in vacuum, and fi is the 
Planck constant divided by 27r, the pair creation rate per 
unit phase space volume given by q{t, r,p) = qo{t, r)F{p), 



where the polarization current is jpoi — EqQ/\E\^ and 
the vacuum polarization charge density is defined via the 
continuity equation dtPpoX + V • jpoi = 0. The particu- 
lar choice of jpoi gives overall energy conservation to the 
system where the energy of newly created pairs is com- 
pensated by a decrease of the electrostatic energy. The 
fluid equation for the cold electrons, positrons, and ions 
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are governed by [Ij, [1 
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valid for length scales larger than vp/^pe, being the 
Fermi velocity and Wpe = (noe^/eoWe)"^^^ is the elec- 
tron plasma frequency for some typical electron density 
no- Here we have added the recombination term for the 
sake of generality, and the index e {p) denotes the elec- 
trons (positrons), i denotes the ions, ne,p,i is the num- 
ber density, Pe,p,i is the fluid momentum, — e (-l-e) is 
the electron (positron) charge, ZiC is the ion charge, and 



(1 + p i/'Tig p iC^Y^'^ is the relativistic gamma 



e,p,i 



factor. Moreover, the generalized relativistic Bohm po- 
tential is defined according to (neglecting the spin of the 
particles) 

UBe,p = ("•e,p/7e,p)"^^^V^(rie,p/7e,p)^/^ 

-(ne,p/7e,pc")-'/'52(ne,p/7e,p)'/^ (3) 

obtained using a Madelung representation for the 
Klein-Gordon equation, and we have defined Pe^p = 
J Pe,pF(fp. Equation l(^ - lf2h|) . together with Eq. 
determines the semi-classical dynamics of the three com- 
ponent plasma where quantum mechanical correction for 
the electrons and positrons are taken into account. 



III. ONE-DIMENSIONAL MODEL 

We now specialize to a one-dimensional geometry. In 
addition to the effects discussed in the 3D model, we will 
assume that the electron-positron pairs are created at 
some distance from each other so that the energy of the 
particles are compensated by a decrease of the electric 
field energy. The Poisson equation is on the form 



dxE ^ —{rip - Ue- ZiUi) 



(4) 



where ions are taken to be immobile so that ZiUi = Uq. 
We assume that the the electron and positron continuity 
equations, respectively, are on the form 



and 
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and the momentum equations 

dtPe + nieC^dxJe = 



-eE 



and 



dtPp + mpC^dx^p = eE, 
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where we have assumed that the cold particles are created 
with momenta equal to the fluid momenta of the electrons 
and positrons, respectively, so that the sources do not 
contribute to the momenta of the particles but only to 
the particle number densities. In this manner, we do not 
have heating of the particles and do not need to include 
pressure or take into account kinetic effects. The source 
terms are taken to be on the form 



= qo + dx { Xjeqo 
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and correspond roughly to that the creation of electrons 
and positrons at the position x depends on the electric 
field at a distance X-jpEcrit/ E and —XjeEcrit/E, respec- 
tively, from X. For a stronger electric field, the the dis- 
tance of creation decreases, while it increases for larger 
7j . This form on the source term gives an approximate 
energy conservation law of the form 
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where the right-hand side can be assumed to be small if 
X\dx\ ^ 1, and where we have assumed that all fields and 
velocities vanish at |x| = oo. The conservation law was 
obtained by using Amperes law 
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which is derived by differentiating the Poisson equation 
by time and using the continuity equations for the parti- 
cles. 
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IV. THE NORMALIZED SYSTEM OF 
EQUATIONS 

We may normalize the system of equations ^ by in- 
troducing the Compton wave length A and the time scale 
T = h/nieC^, normalizing the electric field by -Ecrit, nor- 
malizing all momenta by meC (with mp = me), and nor- 
malizing the number densities by some preferred density 
riQ. Assuming stationary ions, we then obtain the nor- 
malized and dimcnsionless system 



dtUe + V • 

dtUp + V 
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and we have used the stationary ion background den- 
sity as the normalization density hq, Nq = noh^/m^c^, 
the electron plasma frequency is normalized by 1 /r such 
that Wpe = {2aNoy/^/2TT, a = e^/Aneohc 1/137 is 
the fine structure constant, and we have denoted a — 
No{a/X^){ve/c). 



NUMERICAL ANALYSIS OF THE 
ONE-DIMENSIONAL SYSTEM 




FIG. 1: The electrostatic field (left panel) electron density 
(middle panel) and positron density (right panel) as a function 
of space X and time t. We used A'^o = 0.2 and the initial 
conditions Ue = 1.01 -I- 2{x/ L) exp{-x^ / L^) with L = 6000, 
and rip = 0.01, for the electron and positron number densities, 
respectively. 
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FIG. 2; The total number of created electron-positron pairs, 
SNp = SNe = Ns — Ne^t^o, where A^e ~ rie dx, as a function 
of time. We see a decrease in the pair creation rate, correlated 
with a decrease of the amplitude of the electrostatic field seen 
in Fig. 1. 



and are solved numerically. The results are displayed in 
Figs. (l)-(3). Initially, the electron density is perturbed 
locally on the form Ue = 1.01 + 2{x/L) exp(— x^/L^) with 
L = 6000, and for the positrons we take a small non-zero 
component Up = 0.01. All particle species are assumed 
initially to be in rest, i.e., Pe = Pp = at t — 0. We 
have assumed a dense plasma so that A^o = 0.2. Ini- 
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The one-dimensional system of equations are 
normalized so that 

dxE ^ujl^{l-ne+np), (15a) 

d,n, + (^^) =90 + 9. (^) , (15b) 

dtUp + doo [ ^^^^ ) = qo-dx ( ) , (15c) 

\ Ip J \ E J FIG. 3: Thetotalenergy £:tot = ^(ne7e-^ rip 7p-^£;V2^^p) rfa; 

as a function of time (solid curve). The dotted curve displays 

)the same integral for the case when the correction for the 
— —E^ (15d) spatial displacement of the created pairs, the last terms in 

the right-hand sides of Eqs. llSbll and (llScll . are neglected. 
\ In the latter case, the creation of electron-positron pairs lead 

H dx j Pp = E. (15e) to an increase of the total energy in the system. 
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tially, the maximum electric field is approximately half 
the Schwinger field, E^ax — 0.5, and we see in the right- 
hand panel of Fig. 1 that electron-positron pairs are 
created so that the particle density of the positrons are 
increased in a small region around x — where the 
electric field amplitude has its maximum. Due to the 
ultra-strong electric field, the electrons are also accel- 
erated so as to neutralize the plasma, and hence the 
electric field strength decreases with time. We see in 
Fig. 2, that the pair creation rate is largest initially, 
when the electric field is strongest, and decreases at later 
times. Both the electrons and positrons are accelerated 
to ultra-relativistic speeds with gamma factors of the or- 
der 100-500, and the dynamics is essentially a balance 
between the kinetic energy of the particles and the po- 
tential energy stored in the electric field. We consider 
the energy balance in Fig. 3, where we have plotted the 
time evolution of the sum of the total relativistic particle 
energies and the electrostatic energy. We see that this 
energy is approximately conserved (solid line) when the 
the full expressions for the pair creation in the continu- 
ity equations (|15b|l - (|15c|l are used. The last terms in 
the continuity equation account approximately for that 
the electron and positron is created at some spatial dif- 
ference from each other so that the total energy of the 
newly created pair is compensated by a decrease of the 
electric field and hence of the electrostatic energy. If the 
last terms in the right-hand sides of the continuity equa- 
tions are neglected, meaning that the pairs are created at 



the same point in space, then there is a visible increase of 
the total energy when the pairs are created (the dotted 
curve in Fig. 3). However, the energy loss to the pairs is 
relatively small compared to the total free energy of the 
system, stored initially in the electrostatic field. 

In summary, we have presented a dynamical and self- 
consistent model for the electron-positron pair creation 
by a strong electrostatic field in a dense plasma. The 
electrostatic field is excited self-consistently by large- 
amplitude electron waves, which gives rise to electric 
fields that are comparable with Schwinger's critical field. 
We have derived an approximately energy-conserving 
one-dimensional model that takes into account the de- 
crease of the electrostatic energy as electron-positron 
pairs are created. The model presented here consti- 
tutes a first step of understanding the dynamics of large- 
amplitude electrostatic waves that arc strong enough to 
create electron-positron pairs. A future model should be 
derived from first principles of quantum electrodynam- 
ics, to take into account the distribution of the created 
electron positron pairs in momentum space and the self- 
consistent energy conservation of the system. 
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